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We continue the study of a local, gauge invariant Yang-Mills action containing a mass parameter, which 
we constructed in a previous paper starting from the nonlocal gauge invariant mass dimension two operator 
F /JV (D 2 )^ 1 F^. We return briefly to the renormalizability of the model, which can be proven to all orders of 
perturbation theory by embedding it in a more general model with a larger symmetry content. We point out 
the existence of a nilpotent BRST symmetry. Although our action contains extra (anti)commuting tensor fields 
and coupling constants, we prove that our model in the limit of vanishing mass is equivalent with ordinary 
massless Yang-Mills theories. The full theory is renormalized explicitly at two loops in the MS scheme and all 
the renormalization group functions are presented. We end with some comments on the potential relevance of 
this gauge model for the issue of a dynamical gluon mass generation. 
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I. INTRODUCTION. 



Yang-Mills gauge theories, with quantum chromodynamics (QCD) modeling the strong interaction between elementary par- 
ticles as one of the key examples, are quite well understood at very high energies. In this energy region, asymptotic freedom 
QJ | CI S HI 01 se t s in, which in turn ensures that the coupling constant g 2 is small enough to make a perturbative expansion in 
C"| , powers of g 2 possible. The elementary QCD excitations are the gluons and quarks. 

Our current understanding of non-Abelian gauge theories is still incomplete in the infrared region. At lower energies, the 
interaction grows stronger, preventing the use of standard perturbation theory to obtain relatively acceptable results. Nonper- 
turbative aspects of the theory come into play. The most notable, yet to be rigourously proven nonperturbative phenomenon, is 
the fact that the elementary gluon and quark excitations no longer belong to the physical spectrum, being confined into colorless 
■ - - i states such as glueballs, mesons and baryons. 

A widely used strategy to parametrize certain nonperturbative effects of the theory amounts to the introduction of so called 
condensates, which are the expectation values of certain operators in the vacuum. Furthermore, one can employ the operator 
product expansion (OPE) (viz. short distance expansion) which can be applied to local operators, in order to relate the associated 
condensates to nonperturbative power corrections which, in turn, give additional information next to the perturbatively calculable 
contributions. 

As we are considering a gauge theory, these condensates should be gauge invariant if they are to enter physical observables. 
This puts rather strong restrictions on the possible condensates, the ones with lowest dimensionality are the dimension three 
quark condensate (\j/V|/) and the dimension four gluon condensate (F 2 V ). There is a variety of methods to obtain estimates of 
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these condensates, such as the phenomological approach based on the S VZ sum rules, for a recent overview, the use of lattice 
methods, |6], as well as the use of instanton calculus, fin . 

Gauge condensates are necessarily nonperturbative in nature, as gauge theories do not contain a mass term in the action due 
to the requirement of gauge invariance. However, through nonperturbative effects, a nontrivial value for e.g. can arise. 

In (8J], it was already argued that also gauge variant condensates could influence gauge variant quantities such as the gluon 
propagator. In particular, the dimension two gluo n condensate (Aj,) has received much attention in the Landau gauge 
El III [I! [ll IB El over the past few years. An OPE 

argument based on lattice simulations has provided evidence that this condensate could account for quadratic power corrections 



of the form ~ ^y, reported in the running of the coupling constant as well as in the gluon propagator [ 1 lU 1 3Ul 4124 1271 131I.I33I 

This nonvanishing condensate (Afy gives rise to a dynamically generated gluon mass, lEll Il2l fl6l Il9l EoTl . The appearance 
of mass parameters in the gluon two point function is a common feature of the expressions employed to fit the numerical data 
obtained from lattice simulations, l35l l3fil l37tl . Let us mention that a gluon mass has been found to be useful also in the 
phenomenological context, j^Eal- 

The local operator A^ in the Landau gauge has witnessed a renewed interest due to the recent works 0,E3l 5 as me quantity 

which is gauge invariant due to the minimization along the gauge orbits, could be physically relevant. In fact, as shown in 
in the case of compact three-dimensional QED, the quantity (A^ n \ seems to be useful in order to detect the presence 
of nontrivial field configurations like monopoles. One should notice that the operator A^ in is highly nonlocal and therefore it 
falls beyond the standard OPE realm that refers to local operators. One can show that A^ in can be written as a infinite series of 

0(A 4 ). (2) 



nonlocal terms, see 1 40 ,41] and references therein, namely 



^min — 2 ' ^ ' 



However, in the Landau gauge, d^A^ — 0, all nonlocal terms of expression (|2} drop out, so that A^ in reduces to the local operator 
A^, hence the interest in the Landau gauge and its dimension two gluon condensate. However, a complication, as already outlined 
in our previous paper |40], is that the explicit determination of the absolute minimum of A^ along its gauge orbit, and moreover 
of its vacuum expectation value, is a very delicate issue intimately related to the problem of the Gribov copies Il42ll . We refer to 
1 40] for some more explanation and the original references concerning this point. 

Nevertheless, some nontrivial results were proven concerning the operator A^. In particular, we mention its multiplicative 
renormalizability to all orders of p erturbation theory, in addition to an interesting and numerically verified relation concerning 
its anomalous dimension fl5l 14311 . An effective potential approach consistent with the renormalization group requirements 
has also been worked out for this operator, giving further evidence of a nonvanishing condensate (A^) ^ 0, which lowers the 
nonperturbative vacuum energy lfl2il . 

A somewhat weak point about the operator A^ in is that it is unclear how to deal with it in gauges other than the Landau gauge. 
Till now, it seems hopeless to prove its renormalizability out of the Landau gauge. In fact, at the classical level, adding ([0 to the 
Yang-Mills action is equivalent to add the so-called Stueckelberg action, which is known to be not renormalizable 11441 14511 . We 
refer, once more, to [40] for details and references. 

In recent years, some progress has also been made in the potential relevance of dimension two condensates beyond the Landau 
gauge. We were able to prove the renormalizability of certain local operators like: A^ in the linear covariant gauges, (jA^A^ + 

ac"c") in the nonlinear Curci-Ferrari gauges, and (jA^A^ + ac^c-P) in the maximal Abelian gauges jg4]. A renormalizable 
effective potential for these operators has been constructed, giving rise to a nontrivial value for the corresponding condensates, 
and a dynamical gluon mass parameter emerged in each of these gauges 1 46, 47, 48, 49, 50]. There also have been attempts to 
include possible effects of Gribov copies fl8l EH l5ll l52il . Unfortunately, the amount of numerical data available from lattice 
simulations is rather scarce in the aforemention ed g auges. Nevertheless, let us mention that a dynamical gluon mass in the 
maximal Abelian gauge has been reported in 1 53, 54]. In the Coulomb gauge too, the relevance of a dimension two condensate 
has been touched upon in the past 1, 5511 . 

Although the renormalizability of the foregoing dimension two operators is a nontrivial and important fact in its own right, 
their lack of gauge invariance is a less welcome feature. Moreover, at present, it is yet an open question whether these operators 
might be related in some way to a gauge invariant gluon mass. 

Many aspects of the dimension two condensates and of the related issue of dynamical mass generation in Yang-Mills theories 
need further understanding. An important step forwards would be a gauge invariant mechanism behind a dynamical mass, 
without giving up the important renormalization aspects of quantum field theory. 

We set a first step in this direction in our previous paper 140(1 . As a local gauge invariant operator of dimension two does not 
exist, and since locality of the action is almost indispensable to prove renormalizability to all orders and to have a consistent 
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calculational framework at hand, we could look for a nonlocal operator that is localizable by introducing an additional set of 
fields. As pointed out in 1 40], this task looks extremely difficult for the operator A^ in if we reckon that an infinite series of 
nonlocal terms is required, as displayed in Q. Instead, we turned our attention to the gauge invariant operator 



f\f = — I d 4 xF" 
D 2 VT J ^ 



ab 

f& , (3) 



where D 2 = D^D^ is the covariant Laplacian, being the adjoint covariant derivative. The operator already appeared in 
relation to gluon mass generation in three-dimensional Yang-Mills theories ll56Tl . 

The usefulness of the operator (0 relies on the fact that, when it is added to the usual Yang-Mills Lagrangian by means of 
— \rn 2 F jpF , the resulting action can be easily cast into a local form by introducing a finite set of auxiliary fields 14011 . Starting 
from that particular localized action, we succeeded in constructing a gauge invariant classical action S c i containing the mass 
parameter m, enjoying renormalizability. This action S c i was identified to be 

Sphys = S c l + Sgf , (4) 



S c i = d x 



(5) 



32 v " v " v/ 16 

S gf = Jd A x{^b%"+b a d,Al+^d,Dfc b ). (6) 



We notice that we had to introduce a new quartic tensor coupling X ahcd , as well as two new mass couplings A-i and X3. The 
renormalizability was proven to all orders in the class of linear covariant gauges, implemented through S g f, via the algebraic 
renormalization formalism 15711 . Without the new couplings, i.e. when A-i = 0, X3 = 0, \ abcd = 0, the previous action would not 
be renormalizable. 

In this paper, we present further results concerning the action (0 obtained in [40] . In section II, we provide a short summary 
of the construction of the model (|5} and we present a detailed discussion of the renormalization of the tensor coupling X" , not 
given in [40]. We draw attention to the existence of an extended version of the usual nilpotent BRST symmetry for the model 
@. We introduce a kind of supersymmetry between the novel fields {B^,,!}^, G^G^,} which is enjoyed by the massless 
version, m = 0, of the action 0. In section III, we discuss the explicit renormalization of several quantities. The fields and the 
mass m are renormalized to two loop order in the MS scheme. The p" 6 "'-function of the tensor coupling \ abcd is determined at 
one loop, by means of which it shall also become clear that radiative corrections (re)introduce anyhow the quartic interaction in 
the novel fields in the action @. These fields are thus more than simple auxiliary fields, which appear at most quadratic ally. A 
few internal checks on the results are included, such as the explicit gauge parameter independence of the anomalous dimension 
of m. It is also found that the original Yang-Mills quantities renormalize identically as when the usual Yang-Mills action would 
have been used. This is indicative of the fact that the massless version of (|4} might be equivalent to Yang-Mills theory, quantized 
in the same gauge. This is a nontrivial statement, due to the presence of the term proportional to the tensor coupling \ abal . In 
section IV, we use the aforementioned supersymmetry to actually prove the equivalence between the massless version of @ 
and Yang-Mills theories. In the concluding section V, we put forward a few suggestions that might be useful for future research 
directions. 

II. SURVEY OF THE CONSTRUCTION OF THE MODEL AND ITS RENORMALIZABILITY. 

In this section we present a short summary of how we came to the construction of our model (0} in ll40ll . We shall also point 
out a few properties of the corresponding action not explicitly mentioned in H4QI1 . 



A. The model at the classical level. 

We start from the Yang-Mills action Sym supplemented with a gauge invariant although nonlocal mass operator 

Sym + So , (7) 

with the usual Yang-Mills action defined by 



Sym=\J d 4 xF^ v , (8) 



4 



and with 



So 



m 
T 



d 4 xF" 



//V 



rib 



The field strength is given by 



and the adjoint covariant derivative by 



( D r l 

F^ = d^-d v A^-gf^A% 



//V 



ah 



. d5 a b _ gf abc A c 



(9) 



(10) 



(11) 



In order to have a consistent calculational framework at the perturbative level, we need a local action. To our knowledge, it is 
unknown how to treat an action like Q, such as proving its renormalizability to all orders of perturbation theory. This is due 
to the presence of the nonlocal term (|9). As we have discussed in 1 40], the action can be localized by introducing a pair 
of complex bosonic antisymmetric tensor fields, (B^ V ,B^ V ), and a pair of complex anticommuting antisymmetric tensor fields, 



(G A , V ,G^ V ), belonging to the adjoint representation, according to which 



DBDBDGDG exp 



d A x%,D ab D b ^- 



d 4 xG fjV DfD b J c G c fJV ■ 



im 
~4 



fx(B-B) a I% 



(12) 



It is worth mentioning the special limit m = 0, in which case we have in fact introduced nothing more than a rather complicated 
unity written as |85] 



/ 



DBDBDGDGexp 



[- / d«xB^D b a c B^-~ 



1 



d'xG^DfD^ 



= 1 



Hence, we have obtained a local, classical and gauge invariant action 

S = Sym + Sbg + S,i 

where 



(13) 



(14) 



Sbg 

Sm 

The gauge transformations are given by 



1 

4 

im 



d 4 x ( B"D" b Di c B c ;,- 



s-< a T\ab r\bc *~ic 
°//V 1J a u a ^/v 



(15) 





= -Dfa b , 




= gf'^B^ 




= gf' bc co% v 


5G^ V 


= gf bL \o b G^ 




= gf' bc a% v 



with co" parametrizing an arbitrary infinitesimal gauge transformation, so that 

85 = 8 (Sym + Sbg + S m ) = . 



(16) 



(17) 



B. The model at the quantum level. 

Evidently, the construction of the classical action dl4> is only a first step. We still need to investigate if this action can be 
renormalized when the quantum corrections are included. This highly nontrivial task was treated at length in 1 40] to which we 
refer the interested reader for background information. Nevertheless, we shall take some time to explain the main idea as well 
as to present a detailed analysis of the quantum corrections affecting the quartic tensor coupling X cd . 

As the quantization of a locally invariant gauge model requires the fixing of the gauge freedom, we shall employ the linear 
covariant gauge fixing from now on, as it was done in 1 40], and which is imposed via 0. 
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To actually discuss the renormalizability of (I14L we found it useful to embed it into a more general class of models described 
by the action 



Sym + S gf + J d 4 x (TiDfD^Bl - ^DfD b ;Gi 



d\ ( (U^Gf + V'mvB- - Vi^Bf + U ipV G") + (V^ v d 2 V^ v - U ipV d 2 U ifJV ) ) 



d A X%2 (Vi/jvd/jdaViva - ^v^a^/va) - J d 4 xC, (U i M vU i/jV U ja pUj a p 
+ V^ v Vi pV Vj a pVj a p-2Uj MV U ilJV Vj a pVj a p) , (18) 

where use has been made of a composite Lorentz index i = (ju,v), i = 1 ...6, corresponding to a global U (6) symmetry [40] of 
the action Jl 81 . The quantities V^v, Vipv, f/^v and t/ !lUV are local sources. The identification between objects carrying indices i 
and (jU,v) is determined through 



( a" n a r"> n a \ ^ (n a n a r"> n" \ 



(VifiV:Vi[.iViUi l j\,,Ujiiv^ — — (Vopjuvj ^ap^Vi t^ap^Vj t^op/jv) • (19) 

The free parameters Xi» X2 an d C are needed for renormalizability purposes. As far as we are considering Green functions of 
elementary fields, their role is irrelevant as they multiply terms which are polynomial in the external sources. 

The reason for introducing the action (1 1 81 is that for m = 0, the action H41 enjoys a few global symmetries which are lost for 
m ^ 0, whereas the action 1)18} also enjoys these symmetries when the global symmetry transformations are suitably extended 
to the sources. We refer to 1 40] for the details. Evidently, the general action Jl 8I > must possess the action MAI we are interested 
in as a special case. The reader can check that the connection is made by considering the "physical" limit 

\\mV mN = limVo PjU v = (So/iSpv-SovSp^,) , 

phys phys 1 

\\mUappv = limf/op^v = , (20) 

phys phys 

i.e. 

limE = S. (21) 

phys 

In |40j], it was shown that the action dl 81 obeys a large set of Ward identities. We shall not list them here, but mention that the 
action E is invariant under a nilpotent BRST transformation s, acting on the fields as 





= -D ab c b 




sc a 


= if bc c b c c : 




*B% 


= gf abc c b B M ^ 


-G%, 




= gf be c% v , 






= gf^ttCfa , 




sGfjy 


= gf abc c b G^- 




st' 


= b a , 




sb" 


= o, 





(22) 



and on the sources as 









= o, 








= o, 



(23) 
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such that 



sL = 0. 
s 2 = 



(24) 



By employing the algebraic renormalization technique 15711 . it was found that Jl 81 is not yet the most general action compatible 
with all the Ward identities, including the Slavnov-Taylor identity associated to the BRST invariance described by J24i . The 
most general and renormalizable action was identified to be 



(25) 



where 



U (b"b" —~g"g") (Vj^V jtn - U jlN U 



\ abed 



16 



B*B b -~3*G b ) (tfj-Bj -WjGj 



A-3 (B; GjVjjjvU jfj V + G, G"UiuvU i u v + B, B"Vj uV V , uV — G, B"U, u vV ; 



G"B"U ijj V V jpv + Gj BjUifwVjfM ■ 



-B"B a jVj pV V jflV 



-G"G"jU iiJV U^ v 



-B[ BjVj/jvVj/jv + -Gj GjUifiyU j/j V 



(26) 



The quantities X\ and are independent scalar couplings, whereas K a is an invariant rank 4 tensor coupling, obeying the 
generalized Jacobi identity 



^manymbed ^mbnyamcd 



I jrmcny^cibmd j ^mdnyabcm 



0, 



(27) 



and subject to the following symmetry constraints 



yabed ycdab 

yabed ybacd 



(28) 



which can be read off from the vertex that \ abcd multiplies. When we specify the action (125 \ to the physical values (I20> . we 
obtain the main outcome of the paper 14011 . which is the physical action S p h ys given in @, that is renormalizable to all orders of 
perturbation theory in the linear covariant gauge, imposed via S g f. The renormalizability is of course ensured as @ is a special 
case of the more general renormalizable action d25i . since 



limS 

phys 



gen 



Sphys 



(29) 



We notice that the couplings A-i and X3 are in fact part of the mass matrix of the fields {Z>^ V ,Z>^ V ,G^ V ,G^ V }. 

We end this subsection by mentioning that the classical action S c / is also invariant with respect to the gauge transformations 
dl6> . since the terms °= {\\,\3,/\. abcd } are separately gauge invariant. 

C. The renormalization of the tensor coupling X abcd . 

As already mentioned, this section is devoted to providing further details of the renormalization of the tensor coupling X abcd , 
not fully covered in j40ll . 

The term we are interested in at the level of the bare 18611 action is given by 



tt'x 



' T, ab ed 



o / n a nb y-fi y-^b \ ( n c nd f c fd 

\ tS i,o ts i,o - u i,o U i,o J \ tS i,o ts j,o - ij j,o U j,o 



(30) 



where, in the notation of I4 1 



{B,B,G,G} a o . = ^%{B,B,G,G}" = 



1 +T| [ 03 + -«4 



{B,B,G,G}; 



(31) 



where 03 , 04 are arbitrary coefficients and r| stands for a perturbative expansion parameter 1 40] . 
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The most general counterterm corresponding to the renormalization of the 4-point vertex {B"B b — Gf G\ ) (jljB'j — G C jG d 
turns out to be given by 



abcd 



(4« 3 + a 6 ) (b?2J? - 5? Gf ) (l^ - G*G?) , (32) 

where as is a free coefficient and fW abcd is an arbitrary invariant rank 4 tensor, composed of all the other available tensors (such 
as X" bcd , 5 afo , / and invariant objects constructed from these and T"). By the Ward identities, it is nevertheless restricted by 

cy[ abed _ ^ cdab ^ 

M abcd = M bacd , (33) 

which are of course the same symmetry constraints as those for X a , see d28i . Also the Jacobi identity 1271 applies to M abcd . 
The counterterm is thus not necessarily directly proportional to the original tensor X abcd '. This has a simple diagrammatical 

explanation, as diagrams contributing to the 4-point interaction « (B"B b — G"g'-) yEjB'j — G C jG d ^j can be constructed with the 

other available interactions. This also means that, even if X d = 0, radiative corrections shall reintroduce this 4-point interaction. 
This shall become more clear in section III, where the explicit results are discussed. 
We can thus decompose the bare tensor coupling X abcd as 

jabed =z yabcd +z abcd ^ (34) 

where Z and Z abcd contain the counterterm information, more precisely Z contains the counterterm information directly propor- 
tional to X abcd , while Z abcd contains, so to say, all the rest. Evidently, Z abcd will obey analogous constraints as given in d28l > or 
fl33t . The tensor M abcd can be decomposed similarly into 

abcd _ y^abed _|_ ^ abcd _ fcibcd ^5) 

=9\£ abcd 

In the previous paper 1 40], we erroneously omitted the <^[ abcd part. Using ( I32l i. ( I34l i and ( I35l l allows for a simple identification, 
being 

Z = 1 +T|(fl 6 — 2fl4) , 

z abcd = y} (4a 3 + a 6 )K abcd . (36) 

Consequently, the model is still renormalizable to all orders, although X abcd is not multiplicatively renormalizable in the naive 
sense. The situation can be directly compared with Higgs inspired models like the Coleman- Weinberg action l60H6lll . in the 
sense that also there a similar mixing occurs between the different couplings, in casu the gauge coupling e 2 and the Higgs 
coupling X. This is nicely reflected in the p-functions for the couplings, which are series in both e 2 and X. It is even so that 
setting the Higgs coupling A, = does not make p^(e 2 ,A.) vanish. See 16111 for the three loop expressions. As we shall discuss 
later in this paper, the p" &(;rf -function of the tensor coupling X abcd will be influenced by the gauge coupling g 2 . Vice versa, one 
might expect that X abcd could enter, in a suitable colorless combination, the p g 2-function for g 2 . This is however not the case. 
We shall present the general argument behind this in section IV. The p g 2 -function remains thus identical to the well known 
P(g 2 )-function of Yang-Mills theory. 

Let us end this subsection by mentioning that the method of using the extended action (1251 . which is a generalization of 
another action like @ and which exhibits a larger set of Ward identities, turns out to be a powerful tool in order to establish 
renormalizability to all orders. This is reminiscent of Zwanziger's approach to prove the renormalizability of a local action 
describing the restriction to the first Gribov horizon l58tl59ll . 



D. A few words about the BRST symmetry and a kind of supersymmetry. 

Let us now return for a moment to the action S p i iys given in @. As it is a gauge fixed action, we expect that it should have 
a nilpotent BRST symmetry. However, one shall easily recognize that the BRST transformation s as defined in J22i no longer 
constitutes a symmetry of the action (0}. This is due to the fact that setting the sources to their physical values i20\ breaks the 
BRST s as the transformations fl23i are incompatible with the desired physical values ( 12m . 
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Let us take a closer look at the breaking of this BRST transformation s. Let us define another transformation s at the level of 
the fields by 







Jc a 


Q i i 

5 rabc a b 

= V ' 


TFt" 


- gf abc r b R c 




= gf bc c% v 




= gf bc c»G% 


sG^ v 


= gf bc c>% v 


Ic 1 ' 


= b", 




= 0. 



A little algebra yields 

s'Sphys = , (38) 
J 2 = 0. (39) 

Hence, the action S p i iys is invariant with respect to a nilpotent BRST transformation J. We obtained thus a gauge field theory, 
described by the action S p h ys , 0, containing a mass term, and which has the property of being renormalizable, while nevertheless 
a nilpotent BRST transformation expressing the gauge invariance after gauge fixing exists simultaneously. It is clear that Fstands 
for the usual BRST transformation, well known from literature, on the original Yang-Mills fields, whereas the gauge fixing part 
S g f given in can be written as a J-variation, ensuring that the gauge invariant physical operators shall not depend on the 
choice of the gauge parameter |57|. 

We can relate J and s. Let us start from the original localized action (114-1 and let us set m = 0. Then it enjoys a nilpotent 
"supersymmetry" between the auxiliary tensor fields {B^ V ,Z?", V ,G^ V ,G A , V }, more precisely if we define the (anticommuting) 
transformation 8 S as 





— c a 

> 


§.vG^ v 


= 0, 


§sG A , v 


r" 




= o, 



8^ = for all other fields ¥ , (40) 



then one can check that 



*Mm=o) = °- (41) 

Let us mention for further use that, 8 S being a nilpotent operator, it possesses its own cohomology, which is easily identified with 
polynomials in the original Yang-Mills fields {A^,b",c a ,c"}. The auxiliary tensor fields, {B^ v ,B" lV ,G^ v , G|, v }, do not belong to 
the cohomology of 8 S , as a consequence of the fact that they form pairs of doublets 115 711 . 
Taking a closer look upon eqns. J22i . d37i and (I40> . one immediately verifies that 

s = 7+8 s . (42) 

When m ^ 0, the action (I14> is no longer 8 s -invariant. Nevertheless, this 8. s -symmetry can be kept if the more general action 
(1181 is employed, when we extend the 8 4 -invariance to the introduced sources |87] as 









= o, 


Osl/ i/jv 


= v, MV 




= 0. 



Eventually, the most general and renormalizable action J25l > turns out to be compatible with the 8 s -invariance too, as it should 
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be. This is obvious if we recognize that we can write 

S gen = SYM + S.f + bs j d A xG a i Dl b D'^B c i 



+ b s Jd 4 x hiBfG'KVj^Vj.v-Uj^Uj^)) 

+ X abcd ((W -GfGf) (G>$)) +X 3 (BfjU jlM {By^-G-U^)) 

- h 3 (BfU Jllt {GfUi^ + BfVfr)) + hi (G"Vj pV (G-U ifjv - B°V m v)) 



(44) 



and invoke the nilpotency of 8 S . 

What happens when we return to the physical action @? Clearly, the 8 4 -invariance is broken as d20l > and d43t are incompatible. 
The presence of the mass m thus breaks the supersymmetry 8 S . As a consequence, by keeping J42i in mind for the fields, the 
BRST tranformation s is lost too. Fortunately, we recover another BRST invariance 5, d37i . for the physical action. We shall 
come back to the relevance and use of the 8 s -supersymmetry in section IV. 



E. Intermediate conclusion. 



The classical gauge invariant action S c i can be quantized in the linear covariant gauges, whereby a nilpotent BRST symmetry 
and renormalizability to all orders of perturbation theory are present. 

The most famous gauge models exhibiting renormalizability with the possibility of massiv e gauge bosons are of course those 
based on the Higgs mechanism, which is related to a spontaneous gauge symmetry breaking 1 64, '65U661I671 16811. 

Few other Yang-Mills models exhibiting mass terms for the gauge bosons exist. We mention those based on the Stueckelberg 
formalism, which give rise to a nonpolynomial action in the extra Stueckelberg fields. However, these models lack renormal- 
izability 1441 145I1 . Other models are based on the works |62, 63] by Curci and Ferrari. Although the resulting models are 
renormalizable, they do not have a classical gauge invariant counterpart, since the mass terms that are allowed/needed by renor- 
malizability are not gauge invariant terms. Typically, the mass term is of the form jm 2 (A^A^ + occ" c") , where a is the gauge 
parameter |88]. Next to the Curci-Ferrari gauges, the special case of the Landau gauge, corresponding to taking the limit a = 
for the Curci-Ferrari gauge parameter, and the maximal Abelian gauges can also be used to built up such renormalizable massive 
models. Unfortunately, these models have the problem of being not unitary a fact related to the lack of a nilpotent 

BRST transformation 1831 . Nevertheless, in the past few years a lot of interest arose in these dimension two operators from 
the viewpoint of massless Yang-Mills theories quantized in a specifi c ga uge. As these operators turn out to be renormalizable 
to all orders of perturbation theory in the specific gauge chosen fj^l47ll48l.l5Cll . a consistent framework can be constructed to 
investigate the condensation of these renormalizable albeit non gauge invariant operators fTplTfillTTl iiilS. This 
has resulted in a dynamical mass generating mechanism in gauge fixed Yang-Mills theories 



III. TWO LOOP CALCULATIONS. 



We now detail the actual computation of the two loop anomalous dimension of the fields and the one loop p-function of the 
tensor coupling \ abcd . In order to deduce the renormalization group functions, there are two possible ways to proceed. One is 
to regard the extra gluon mass operator as part of the free Lagrangian and work with completely massive gluon and localizing 
ghosts throughout. It transpires that this would be extremely tedious for various reasons. First, although the propagators will 
be massive there will be a 2-point mixing between the gluon and localizing ghosts leading to a mixed propagator. Whilst it is 
possible to handle such a situation, as has recently been achieved in a similar localization in fTTll . it requires a significantly large 
number of Feynman diagrams to perform the full renormalization. Moreover, one needs to develop an algorithm to systematically 
integrate massive Feynman diagrams where the masses are in principle all divergent. Although algorithms have been developed 
for similar but simpler renormalizations, we do not pursue this avenue here mainly because the extra effort for this route is not 
necessary given that there is a simpler alternative. This is to regard the mass operator as an insertion and split the Lagrangian into 
a free piece involving massless fields with the remainder being transported to the interaction Lagrangian. Hence to renormalize 
the operator will involve its insertion into a massless Green function, after the fields and couplings have been renormalized in 
the massless Lagrangian. This is possible since it has been demonstrated that the ultraviolet structure of the renormalization 
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constants remain unchanged in MS whether the gluon mass operator is present or not |40]. Moreover, given that the massless 
field approach is simpler and more attractive, we can use the MINCER algorithm to perform the actual computations. This 
algorithm, 1 72], written in the symbolic manipulation language FORM, VI3. FmIi . is devised to extract the divergences from 
massless 2-point functions. Therefore, it is ideally suited to deduce the anomalous dimensions of the fields. Hence we note that 
for the computations the propagators of the massless fields in an arbitrary linear co variant gauge are, |40], 



{c"(p)c b (-p)) 



~n2 



(1-0C) 



PpPv 



5 ab 

~n2 



(G a MG b ap (-p)) 



If- 
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[S,ua8vp — 
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-P)) 



§up §va] 



§up §va] 



(45) 



where p is the momentum. Using Qgraf, l75ll . to generate the two loop Feynman diagrams we have first checked that the same 
two loop anomalous dimensions emerge for the gluon, Faddeev-Popov ghost and quarks in an arbitrary linear covariant gauge as 
when the extra localizing ghosts are absent. This is primarily due to the fact that the 2-point functions of the fields do not involve 
any extra ghosts except within diagrams. Then they appear with equal and opposite signs due to the anticommutativity of the 
G^ v ghosts and hence cancel. This observation shall be given an explicit proof in section IV. We thus note that the expressions 
obtained for the renormalization group functions are the same as the two loop MS results of B [z3 [Hi • For the localizing 
ghosts there is the added feature that the properties of the X abcd couplings have to be used, as specified in J27b and d28i . We have 
implemented these properties in a Form module. However, we note that in the renormalization of both localizing ghosts, we 
have assumed that 



1 



1 
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N A N a 

which follows from the fact that there is only one rank two invariant tensor in a classical Lie group. If this is not satisfied then 
one would require a 2-point counterterm involving the X abcd couplings which was not evident in the algebraic renormalization 
technology which established the renormalizability of the localized operator. Hence, at two loops in MS we find that 

v2 



jB(a,X) = Jg(o,X) = (oc-3)fl + 



— + 2a- — 
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L ,4 



10 
-T F N f 
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i abcdi acbd 



128% 



(47) 



where Na is the dimension of the adjoint representation of the colour group, a = g 2 /(16% 2 ) and we have also absorbed a factor 
of 1/(4tc) into X abcd here and in later anomalous dimensions. These anomalous dimensions are consistent with the general 
observation that these fields must have the same renormalization constants, in full agreement with the output of the Ward 
identities l40ll . 

In order to verify that (147 \ is in fact correct, we have renormalized the 3-point gluon B^ v vertex. Since the coupling constant 
renormalization is unaffected by the extra localizing ghosts (and we have checked this explicitly by renormalizing the gluon quark 
vertex), then we can check that the same gauge parameter independent coupling constant renormalization constant emerges from 
gluon B" N vertex. Computing the 7 one loop and 166 two loop Feynman diagrams it is reassuring to record that the vertex is 
finite with the already determined two loop MS field and coupling constant renormalization constants. Prior to considering the 
operator itself, we need to determine the one loop p-function for the X a couplings. As this is present in a quartic interaction 
it means that to deduce its renormalization constant, we need to consider a 4-point function. However, in such a situation 
the Mincer algorithm is not applicable since two external momenta have to be nullified and this will lead to spurious infrared 
infinities which could potentially corrupt the renormalization constant. Therefore, for this renormalization only, we have resorted 
to using a temporary mass regularization introduced into the computation using the algorithm of 1 78] and implemented in Form. 
Consequently, we find the gauge parameter independent renormalization 



aba! 



= X 



abed 
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(48) 



from both the ^ a/ ^ v B / ^ v Z? a pfi' /op and X clbcd B", v B b ^Gl p G' la P vertices where df cd is the totally symmetric rank four tensor 
defined by 



jabed t> / rj^a^rib ^c^p^ 

a A — LT y 1 A 1 A l A l A 



(49) 
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with T£ denoting the group generator in the adjoint representation, |79]. Dimensional regularization in d = 4 — 2e dimensions is 
used throughout this paper. Producing the same expression for both these 4-point functions, aside from the gauge independence, 
is a strong check on their correctness as well as the correct implementation of the group theory. Unlike for the gauge coupling 
and its p-function, the \ abcd p-function contains terms also involving the gauge coupling g 2 at one loop. Hence, to one loop 
$f cd (a,X) is given by 



- ^x abpq X cpdq + \ a P bc i\ cd P c i + x apcq X bpdq + \"P dc i\ b i ,cc i 

- \2C A \ abcd a + 8C A f abp f cdp a 2 + l6C A f adp f bcp a 2 + 96df cd a 2 ] , 



(50) 



such that in d dimensions 



,, 5 yabcd 
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abed 
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(51) 



Another useful check is the observation that fi abcd enjoys the same symmetry properties as the tensor \ abcc \ summarized in J28t . 

It is worth noticing that X" bcd = is not a fixed point due to the extra A. afc "'-independent terms. Put another way, if we had not 
included the X a/, "'-interaction term in the original Lagrangian, then such a term would be generated at one loop through quantum 
corrections, meaning that in this case there would have been a breakdown of renormalizability. Further, we note that with the 
presence of the extra couplings, the two loop term of this P-function is actually scheme dependent. 

Finally, we turn to the two loop renormalization of the localized operator itself, or equivalently of the mass m. The operator 
can be read off from dl5> and is given by 



(52) 



To do this we extend the one loop calculation, Ecft - by again inserting this operator into a A^-B^ 2-point function and deducing 
the appropriate renormalization constant Z , defined by 



Z n O 



(53) 



One significant advantage of the massless field approach is that there is no mixing of this dimension three operator into the 
various lower dimension two operators, which was evident in the algebraic renormalization analysis, and would complicate this 
aspect of the two loop renormalization. In other words following the path of using massive propagators would have required us 
to address this mixing issue. Hence, from the 5 one loop and 131 two loop Feynman diagrams, we find the MS renormalization 
constant 
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and therefore, 
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(55) 



as the two loop MS anomalous dimension, which is defined as |4< 



y (a,\)=/j—\nZ 



(56) 



As at one loop it is independent of the gauge parameter, as expected from the fact that the operator is gauge invariant. Also, 
the two loop correction depends on the \ abcd couplings as well as the gauge coupling, as expected from our earlier arguments. 
We end this section by mentioning that a factor of (—2) was erroneously omitted in the one loop anomalous dimension y (a) in 
eq.(6.9)of 
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IV. EQUIVALENCE BETWEEN THE MASSLESS THEORY AND USUAL YANG-MILLS THEORY. 

In this section, we shall discuss the usefulness of the 8 s -supersymmetry, defined by (I40i . and show that Green functions which 
are built from the original Yang-Mills fields [90] {A^,c a ,c",b a } are independent from the tensor coupling X" bcd , in the massless 
version of the physical model 

Next to this result, we shall also prove the stronger result that any Green function constructed from the original Yang-Mills 
fields {A^,c a ,c a ,£>"}, evaluated with respect to the massless version of our action, gives the same result as if the Green function 
would be evaluated with the original Yang-Mills action. 



A. The massless case m = 0. 

As we have already noticed in H3\ . the case corresponded originally to the introduction of a unity into the usual Yang-Mills 
action. Evidently, we expect that the model obtained with m = would be exactly equivalent to ordinary Yang-Mills theory. 
Nevertheless, this statement is a little less clear if we take a look at the massless action 



nm=0 
phys 



SyM + Sgf 
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d x 



yjibcd 
~l6~ 
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£> po £) pa "pa^pa 



(57) 



which is obtained from @. The reader shall notice that the quartic interaction °< X abcd is anyhow generated, making the path 
integration over the tensor fields no longer an exactly calculable Gaussian integral. Hence, we could worry about the fact the 
tensor coupling X" bcd might enter the expressions for the Yang-Mills Green functions, which are those built out of the fields A^, 
c a , c a and b", when the partition function corresponding to the action (157ft would be used. 

We recall here that the action (1571 is invariant under the supersymmetry ( 140b . This has its consequences for the Green func- 
tions. Let us explore this now. Firstly, we consider a generic n-point function built up only from the original fields { A" n c" ,c",b"}. 
More precisely, we set 



§ n (x\, . . . ,Xn) — {Gn{Xl, . . . ,Xn))<jm=0 
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DQG n (xi,...,x n )e M»» , 



(58) 



with 



G„{x u . ..,x n ) = A(xi) .. .A(xj)c(x i+ i) . ..c(xj)c(x j+ i). . .c(x k )b(x k+ i) . ..b{x„) 



(59) 



and we introduced the shorthand notation <t> denoting all the fields. We notice that 8 4 G„ = whereas G„ ^ 8 V (. . .), i.e. any 
functional of the form d59l belongs to the 8. s -cohomology. 

We are interested in the dependence of Q n on X abcd . A small computation leads to 
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(60) 



The last line of ibOt is based on the fact that 8 S annihilates the vacuum as it generates a symmetry of the model. 

We have thus shown that all original Yang-Mills Green functions will be independent of the tensor coupling \ abcd . These are 
of course the most interesting Green functions, gauge variant (like the gluon propagator) or invariant (the physically relevant 
Green functions). The previous result does not mean that we can simply set \ abcd = and completely forget about the quartic 



interaction °< (B,,,,B 



(i\ *'( \ w (iv w fiv) (^pa^pa — GpaC?pa J . There is a slight complication, as quantum corrections reintroduce the 



b _n a 



G,vG, 



quartic interaction X abcd even when we set X abcd = 1 9 1 ] . 
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B. Exact equivalence between the massless action and the Yang-Mills action. 

In this subsection, we shall use once more the nilpotent 8 S symmetry to actually prove that 

(G„(xi,...,x„)) SYM+Sgf = {G n (xi,...,x n )) sm =o , (61) 

meaning that the expectation value of any Yang-Mills Green function, constructed from the fields \A",c a ,(?,b a } and calculated 
with the original (gauge fixed) Yang-Mills action Sym + Sgf, is identical to the one calculated with the massless action S"^, 
where it is tacitly assumed that the gauge freedom of both actions has been fixed by the same gauge fixing. 

In this context, let us also mention that the physical content of the massless theory described by SF^Z will not depend on the 

extra tensor fields, as (B'L, G^ v ) and (B^ v , G^ v ) are both 8 s -doublets, and hence any physical operator shall certainly not depend 
on these fields. Physical operators belong to the 8 4 -cohomology, which is independent of 8^ -doublets |57]. More precisely, 
Yang-Mills theory and the massless action ( I57> will have identical physical operators, which belong to the BRST s-cohomology. 
Next to this, we also know that the extra fields of the massless model shall decouple from the physical spectrum as they belong 
to the trivial part of the additional 8 v -cohomology. 

Let us now prove the statement iblt . We shall first prove the following. 

Theorem I Let So be an action constructed from a set of fields (j),, that enjoys a symmetry generated by a nilpotent operator 8. Consider 
a second action Si = So + AS, whereby AS is constructed from the fields (j),, and an extra set (cp^., 7p k ) whereby cp^ and 7p k 
are 8-doublets, such that Si also enjoys the symmetry generated by 8. We assume that the renormalizability of So and Si 
has been established. 

The physical operators of Sq and Si both belong to the 8-cohomology, whereas the extra fields (9^,9^.) do not. This is due 
to the fact that these fields give rise to pairs of 8-doublets, thus having vanishing cohomology. The difference AS = Si — So 
is then also necessarily 8-exact. Let H be an operator belonging to the 8-cohomology. Then we can write 

(#) Si = J d$idy k M e- s °- hs . (62) 

We define a new action 

S k = S + kAS, (63) 



where K is a global 1 92] parameter put in front of the action-part AS. If we multiply the counterterm part, corresponding to 
AS, which belongs to a trivial part of the 8-cohomology, by K, then the renormalizability will be maintained, without the 
need of introducing a counterterm for K, so that Ko = K. The parameter K can be used to switch on/off the difference AS. 
More precisely, we can interpolate continuously between So and Si . Using this, it is not difficult to show that 

due to the 8-exactness of AS and 8-closedness of H . As a consequence 



= (65) 



A related theorem is the following. 



Theorem II Let So be an action constructed from a set of fields (j),. Consider a second action Si = So + AS, whereby AS is constructed 
from the fields (j),, and an extra set (9*,%) whereby 9^ and 9^ are 8-doublets, such that Si enjoys the symmetry generated 
by the nilpotent operator 8. We trivially extend the action of 8 on the fields (j), as 8(|), = 0, so that of course bSo = 0. We 
assume that the renormalizability of So and Si has been established. 

The difference AS = Si — So is necessarily 8-exact. Moreover, the physical operators of Si must belong to the 8- 
cohomology, which is independent of the 8-doublets (9^,9^)- The operators constructed from the fields (]),, i.e. all 
operators of the model So, certainly belong to the 8-cohomology of Si. Let %_ be such an operator. A completely similar 
argument as used in Theorem I allows to conclude that 



Ws = W Sl 



(66) 
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Let us now comment on the usefulness of the previous theorems. Theorem II is applicable to the Yang-Mills action Sym + S g f and 
the massless action S"^, where 8 S is the nilpotent symmetry generator of S™jl®, trivially acting on Sym + S g f. Said otherwise, 
we have just proven that the massless physical model d57l is equivalent with Yang-Mills, in the sense that the Green functions 
of the original Yang-Mills theory remain unchanged when evaluated with the action J57> . 

An important corollary of the previous result is that the running of the gauge coupling g 2 will be dictated by the usual (3(g 2 )- 
function known from common literature, as the renormalization factor for g 2 can be extracted from original Yang-Mills n-point 
Green functions. This result was confirmed in section III, as well as the fact that the other Yang-Mills renormalization group 
functions remain unaltered, again in agreement with Theorem II. 

C. The massive case m ^ 0. 

It might be clear that the 8 s -supersymmetry was the key tool to prove the X" bcd independence in the massless case, as well as 
the equivalence with Yang-Mills theory. 

Evidently, we are more interested in the case that the mass m is present. The question arises what we may say in this case, as 
the supersymmetry is now explicitly broken, i.e. 

8 s S pk ys¥= 0, (67) 

with the massive physical action given in @. The role of the tensor fields {B^ V1 B^ V ,G^ V , G", v } remains an open question, as they 
do not longer constitute a pair of doublets. But, we repeat, gauge or BRST invariance is kept, at both the classical and quantum 
level. 



V. CONCLUSION AND OUTLOOK. 

Returning to the rationale behind the construction of the action @, we recall that is was based on the localization procedure, 
given in J12l> . of the nonlocal operator coupled to the Yang-Mills action, as displayed in @. Clearly, we can return from the local 
physical action to the original nonlocal one only if the extra couplings X\, \$ and X" bcd would be zero. As already explained, even 
if we set these couplings equal to zero from the beginning, quantum corrections will reintroduce their corresponding interactions. 

Moreover, as the massive physical action is not S s -supersymmetric, we cannot simply say that we can choose the extra 
couplings freely. Trying to get as close as possible to the original localized version d!5l > of the nonlocal action @ where only 
one coupling constant is present, we can imagine taking the tensor coupling \" bcd to be of the form 

yabcd = e abcd g 2 + Ma bcd g 4 + ^ (6g) 

i.e. we make it a series in the gauge coupling g 2 , where the coefficients l" bcd are constant color tensors with the appropriate 
symmetry properties as J27b and d28l >. This already avoids the introduction of an independent coupling. We have temporarily 
reintroduced the Planck constant h to make clear the distinction between classical and quantum effects. At the classical level, 
we could try to set ^ cd = in order to kill the quartic interaction. Doing so, we would at least keep the classical equivalence 
between the local and nonlocal actions and @ by employing the classical equations of motion. Unfortunately, this is not 
possible. We should assure that d68i is consistent with the quantum model, meaning that we should assure the consistency with 
the renormalization group equations. Taking the quantum effects into account, next to the divergent contributions canceled by 
the available counterterms, there will be also quantum corrections to the quartic interaction which are finite but nonvanishing in 
Hg 2 . We can fix the classical (tree level) value i c ^ cd by demanding that the proposal J68t i is consistent with the renormalization 
group function ( 15 1> . and likewise for the higher order coefficients t" bcd , i.e. we should solve 

^X ab "<{g 2 ) = p fe 2)^L^V) = W bcd {g 2 ) , (69) 

order by order, with X" bcd (g 2 ) defined as in d68i . By using the renormalization group function d50i it is apparent that £ C Q bcd — 
is not a solution. 

It is interesting to notice that the classical action should already contain the classical quartic coupling X" bcd — l^ cd g 2 in order 
to allow for a consistent extension of the model at the quantum level. This classical value is in fact dictated by one loop quantum 
effects, as it is clear from the lowest order term of d69i . 

A completely similar approach could be used for the mass couplings X\ and A3. We could eliminate the extra couplings in 
favour of the gauge coupling, without making any sacrifice with respect to the renormalization group equations. Of course, this 
is a nontrivial task, as it should be checked whether for instance ( I69> possesses meaningful solution(s), as the coefficient tensors 
gabcd shoujd be a j \ eas i realvalued, whereas the uniqueness of the solution might be not evident. 
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The aforementioned procedure is not new, as we became aware of works like 1 80| l8lll and references therein, where the 
reduction of couplings and its use were already studied. 

In this paper, the role of the extra couplings was not the primary motivation. The main purpose of this paper was to establish 
a further study of the gauge model @ itself. We have reported on a few properties. We mentioned the existence of a nilpotent 
BRST symmetry of the physical action @ and we have commented on the fact that the physical action can be embedded 
into a more general action that possesses an interesting supersymmetry amongst the new fields. We briefly returned to the 
renormalizability, in particular on the algebraic renormalization of the extra tensor coupling \ abcd where we drew an analogy 
with, for instance the Higgs inspired Coleman- Weinberg model. The physical action itself only enjoys the supersymmetry in the 
massless limit, in which case we were able to prove the equivalence with ordinary Yang-Mills theory by using an argument based 
on the cohomology of that supersymmetry. We also presented explicit results concerning the renormalization group functions 
of the model evaluated in the MS scheme: the anomalous dimensions of the original Yang-Mills fields and parameters were 
calculated to two loop order and turned out to be identical to the ones calculated in massless Yang-Mills, supplemented with the 
same gauge fixing, in agreement with the general argument that both models are equivalent in the massless limit. The one-loop 
anomalous dimension of the tensor coupling \ uhcd has also been evaluated, giving explicit evidence that \ abcd = is not a fixed 
point of the model, and finally also the two loop anomalous dimensions of the auxiliary tensor fields as well as of the gauge 
invariant operator — (B — B)F, given in d52i . or equivalently of the mass ;«, have been calculated. 

An interesting issue to focus on in the future would be the possible effects on the gluon Green functions arising from the 
massive physical model (0}, to find out whether mass effects might occur in the gluon sector as, for instance, in the gluon 
propagator. 

Returning to the necessity of introducing extra couplings, since we cannot keep the equivalence between actions and @ 
due to these couplings, the relation with the nonlocal gauge invariant operator F-^F has become obscured. However, this is 
not an unexpected feature, as this operator is highly nonlocal, due to the presence of the inverse of the covariant Laplacian. 
Nevertheless, we believe that the final model @ is certainly relevant per se. To some extent, it provides a new example of 
a renormalizable massive model for Yang-Mills theories, which is gauge invariant at the classical level and when quantized it 
enjoys a nilpotent BRST symmetry. 

There are several other remaining questions concerning our model. At the perturbative level for example, it could be inves- 
tigated which (asymptotic) states belong to a physical subspace of the model, and in addition one should find out whether this 
physical subspace can be endowed with a positive norm, which would imply unitarity. Although the resolution of this topic is 
under study, it is worth remarking that the nilpotency of the BRST operator might be useful in this context. 

Of course, at the nonperturbative level, not much can be said at the current time. The model is still asymptotically free, 
implying that at low energies nonperturbative effects, such as confinement, could set in. One could search for indications of 
confinement, similarly as it is done for usual Yang-Mills gauge theories. An example of such an indication is the violation of 
the spectral positivity, see e.g. I2H I82L l83fl . Proving and understanding the possible confinement mechanism in our model is 
probably as difficult as for usual Yang-Mills gauge theories. 

Finally, it would be interesting to find out whether this model might be generated dynamically. A possibility would be to start 
from the massless version of the action 0, which was written down in J57> . After all, it is equivalent with massless Yang-Mills, 
and we can try to investigate whether a nonperturbative dynamically generated term m(B — B)F might emerge, which in turn 
could have influence on the gluon Green functions. 
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